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Abstract 



We study the weak interaction between a pair of well-separated coherent structures in 
possibly non-local lattice differential equations. In particular we prove that if a lattice 
differential equation in one space dimension has asymptotically stable (in the sense 
of Chow, Mallet-Paret and Shen [2]) traveling wave solutions whose profiles approach 
limiting equilibria exponentially fast, then the system admits solutions which are nearly 
the linear superposition of two such traveling waves moving in opposite directions away 
from one another. Moreover, such solutions are themselves asymptotically stable. This 
result is meant to complement analytic or numeric studies into interactions of such 
pulses over finite times which might result in the scenario treated here. Since the 
traveling waves are moving in opposite directions, these solutions are not shift-periodic 
and hence the framework of Chow, Mallet-Paret, and Shen does not apply. We overcome 
this difficulty by embedding the original system in a larger one wherein the linear part 
can be written as a shift-periodic piece plus another piece which, even though it is non- 
autonomous and large, has certain properties which allow us to treat it as if it were a 
small perturbation. 

1 Introduction 

1.1 The system, hypotheses and main results 

This paper concerns weak interactions between coherent objects in lattice differential 
equations. These interactions include pulse-pulse interactions, the gluing of fronts and 
backs to make a wide pulse, front stacking, and the interaction between a pulse and a 
front. 

We study the equation 



where L G C{X , X) annihilates constant functions and G : X ^ X includes nonlinear 
terms which may be nonlocal, see (HO) below. 

We use II ■ II to denote the norm for this space. 



X = LX + G{X) =: F{X) 



X £X 
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We say that a solution X of ([T]) is a traveling wave if it is of the form Xj{t) = (j){j — ct) 
for some continuous function (f> which has finite hmits at ±00, 

4>{—oo) = a, 0(00) = uj. 

In the case that a = cu, (j) is called a pulse; in the case that a lo we call (j) a front. 

We are interested in proving the existence and stability of solutions which are 
roughly the linear superposition of two separated traveling waves which move with 
different speeds, in particular when the waves are separating apart from one another 
as t increases — a situation we call an exit. Thus we assume that ([1]) admits a pair of 
stable traveling wave solutions 4>~ and each of which is either a pulse or a front. 
We denote their wave speeds and asymptotic values with ± subscripts. We assume that 
c_ < c+ and that 0+ is "located" to the right of . Therefore we require w_ = a+. 
Since L annihilates constant sequences, we can take a+ = without loss of generality. 

To see this, let X -.^ X ~ a and note X LX + G{X) where G'(X) G{X + a). 

Before we can state our main theorem, we need to make precise the hypotheses that 
we impose. In what follows, 

Xb:^{Xe£^ I \\X\\h:=sup\il + e'>^)X,\<oo}, 

the space of functions which decay exponentially fast as j (or — depending on the sign 
of b) goes to infinity. 

Standing Assumptions 

• (HO) (continuity of G) G : X ^ X is oi the form 

G(X)„ = giNi{X)n, ■ ■ ■ Nj{X)„) where g e C/j,'(M"-', R") with .g(0) = and 
Ni G C{Xp) for all (3 G [—6, h] and furthermore commutes with the shift. 

• (HI) (existence of traveling waves) There is a 5 > such that the LDE ^ 
admits traveling wave solutions 0^ e Xb and 0^ G X-b with speeds c_ 7^ c+ and 
LO- = a+. We further assume that 0^ G X^b- 

• (H2) (spectral stability of traveling waves) Let $^(t,to) denote the time t map 
for the linear equation Y = {L + G'{(j)^))Y, let S denote the shift on X, 

{Sx)j = Xj-i and let S'-i$±(l/c±, 0). Then one is a simple eigenvalue of 

(with eigenfunction 0'^) and (j{A^) \ {1} is contained in the open unit disc. 
Here the spectrum is computed regarding A^ as an operator on Xzpb- 

Remark 1.1. Typically, the conjugated operator = (1 + e''^)A±[ ^^^^j ] is a small 
perturbation of A± so long as b is chosen sufficiently small. Thus the spectrum of A'^ 
in X coincides with that of A± in Xb- Hence (H2) may he obtained in examples as a 
consequence of the corresponding stability criterion with X replacing Xb . 

Note that in ^9J, (HI) and (H2) are shown to be sufficient to conclude the asymptotic 
stability of the traveling wave. In Section 11.31 we further discuss these hypotheses as 
they relate to a number of different systems of interest. We now state our main theorem. 

Theorem 1.2. // c_ < c^, lo^ = Q!-|_ and (H0)-(H2) are satisfied, then there exists a 
positive constant a such that for each e > there exist positive constants C , 5q, and t* 
such that if 

\\X,mt - 0+(- - r+) - 0-(- - T_)|| < <5 < 5o 
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with 

T+ — T_ > T 

then there are real constants 7^ and 7^ in (— e, e) such that the solution X of ^ with 
initial condition Xinn satisfies 

e'^'WXit) - - c+t ~T+~ 7,+) c-t - T_ - 7*^)11 < C{e---' + VS) (2) 

for all t > 0. 

We can rephrase this theorem in terms of the "exit manifold" : 

A^e«t - r+) + - T_) : r_ - r+ > r*} . 

■Mexit is a smooth two-dimensional submanifold of £°° (see Proposition 3.4 in [9])and 
consists of all linear superpositions of two well-separated traveling waves. It is not an 
invariant manifold for ([T]) but our main theorem implies that is a local attractor for the 
dynamics. That is: 

Corollary 1.3. If disteaa(Xi„it, -Mexit) < So, then the solution X{t) of ([T]) withX{0) = 
Xtnit satisfies: 

disti^{X{t),Mexu) < Ce-'^* 

Remark 1.4. There are numerous results concerning the existence and stability of 
multipulse solutions for reaction diffusion PDE. For instance: |21 liiSt \26^ treat 

existence and stability of multipulse standing solutions; \13\. \14\ \15\ \25\ \30^ deal 
with counter-propagating fronts and pulses in scalar systems using comparison principle; 
U(A \31}j . handle long distance weak interactions between standing pulses; |5l \27\ 
deal with exit or shooting solutions to systems of reaction diffusion equations, and the 
methods used there are most similar to ours. Additionally, multipulse solutions in a 
Hamiltonian lattice have been studied by the first author in \19\. \2Uj. 

The remainder of this paper is organized as follows. In Section 11.21 we outline our 
approach to the proof of Theorem 11.21 In Section 11.31 we discuss some examples of 
Section [2] decomposes the problem into stable and center eigenspaces and we make 
estimates on this decomposition in Section [3l Finally, Section |4] contains the proof of 
Theorem HH 

Acknovifledgements: The authors would like to express their gratitude to the NSF 
for funding this project under grants DMS 0603589 (AH) and DMS 0807738 (JDW). 
Additionally, special thanks is due to Erik Van Vleck for suggesting this problem. 

1.2 General strategy 

We seek a solution of the form Xj{t) = <p^{j — c^t) + (/)^(j — c+<) + w where w goes to 
zero in as i — > 00. Note that for large values of t, the sum (f)^ {j — c-t) + 0+ {j — c+t) 
is close to zero for compact sets of spatial indices j. To that end, we embed ^ into the 
system 

X- = LA" + G(A-) + H-{t){G{X~ +X+) - G{X-) - G{X+)} =: F_(A-, A+) 

X+ = LX+ + G{X+) + H+{t){G{X- + A+) - G(A") - G(A+)} =: i^+(A-, A+). 

(3) 
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Here and iJ+ are localization operators defined as follows. Let h{x) = for a; < 
and h{x) = 1 for a; > denote the usual Heaviside function, and c = {c-+c+) /2. Define 
the operator H^{t) which acts on spaces of sequences by {H^ {t)X)j = h{i — ct)Xj and 
H^{t) = Id — H^{t). At time these operators localize sequences to the right and left 
half-lattices which are centcrerd "halfway between" 0~ and (p'^ . 

Note that if {X-,X+) solves then X = X' + X+ solves ([J). Thus if a solution 
solves ^ with X^ = 0^ + and X^ = 0+ + w+, where w^{t) are 
decaying to zero, then X ~ X~ + X^ is of the form that we seek with w = + w~ . 
The equation ([3]) is a perturbation of two copies of ([T]). However, the coupling terms 
H^[t){G{X- +X+)-G{X-)-G[X+)} arenotsmah, at least when viewed on A". To 
wit, an application of the Mean Value Theorem shows (roughly speaking) that we have: 

\H-{t){G{X- +X+)-G{X-)-G{X+)]\ < GH-{t)\X-\\X+\ 

< GH-{t) (|0-||0+| + \^-\\w+\ + + \w-\\w+\) 

{t) localizes functions to the left half- lattice, where 0+ G is exponentially small. 
Thus H~\(j)'^\ is exponentially small, and we can handle two of the four terms above. 
The term quadratic, and thus also can be made small. However, H~{t)\(j>~\ 

is 0(1) and thus the term is problem. If w"*" is exponentially localized 

to the right half-lattice, then H~ {t)\w^\ will be small just as was H~\(l)'^\. (We make 
this heuristic argument rigorous in Proposition 13.11 below. ) 

Therefore we will require this localization. For the remainder of the paper we regard 
([31) (after a series of non-trivial changes of coordinates) as an evolution equation in the 
phase space 3^ := x At first blush, this may seem to shrink the size of the space 
of initial data we allow for our equation ([1]). However, for any data Xit^) G A" we have 
H^{to)X{to) e X^b- Therefore we set X^{to) := H^{to)X{to) so that initially X^ and 
Ar+ are supported on the left and right half-lattices respectively. Additionally, for any 
6, Afc C A". Thus the study of ^ in y contains the dynamics of llj in X . (Note that 
for t > to, in general we wiU have X-{t) ^ H-{t)X{t) and X+(t) ^ H+(t)X{t).) 

1.3 Examples 

Note that the class of models ^ which satisfy (HO) is quite general. Namely, any 
system of lattice differential equations satisfies (HO) so long as very mild restrictions on 
the nonlinear piece of the nonlocal coupling are satisfied. The class of LDEs to which 
Theorem 1 1 . 2 1 can be applied is much smaller. In many cases (HI) and (H2) are known 
to be false. Examples are furnished by the Hamiltonian lattices which conserve energy 
and hence cannot satisfy (H2). 

The kind of equations that we have in mind are spatial discretizations of possibly 
nonlocal reaction-diffusion equations. These equations are dissipativc, thus (H2) is not 
immediately ruled out. However, establishing (HI) and (H2) is highly nontrivial. To 
demonstrate this, consider the simple scalar equation 

iln = -^{Un+l + Un-l - 2Un) - f {Un) (4) 

which arises as a spatial discretization of the PDE 

ut = u,j,^ - f{u). (5) 
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Here / is the derivative of a double-well potential, e.g. f{u) = u{u — l)(w — a) for some 
a £ (0, 1). Upon substituting the traveling wave ansatz Un{t) = 4>{n~ ct) we obtain the 
mixed type equation 

- - + 1) + <f>{^ - 1) - 20(0 - fim) (6) 

which is ill-posed as a dynamical system on the infinite-dimensional phase space C([— 1, 1], M). 
Comparing with —c(j)' = — /{(I') which arises as a wave profile equation for the PDE 
([U wc can see why the existence and stabilty theory for Q has lagged behind that for 
([S]). Nevertheless both (HI) and (H2) have been established for The existence the- 
ory (HI) can be based upon either topological fixed point theorems |32j or comparison 
principles [16] . Mallet-Paret has developed the Fredholm theory of differential-difference 
operators |24| and built a continuation argument on this theory [23j which establishes 
(HI) for a more general subclass of ^ than ^ under the mild assumptions of finite 
interaction length, spatial homogeneity, and ellipticity (see [53] for details). The as- 
sumptions of finite interaction length and spatial homogeneity have been weakened in 
[2] and [8]. 

With regards to stability theory, it is usually the case that the essential spectrum 
can be easily computed e.g. via Fourier transform. However, the eigenvalue problem 
is of the form ([6]) with an additional spectral parameter. When comparison principles 
are available, this problem is tractable. When comparison principles are not available, 
little is known. 

We should note that comparison principles are typically available in scalar equations 
of reaction diffusion type and can be used to construct a stable monotone front. Note 
also that if (c, 0(0) is a solution of equation ([5]), then (— c, (/)(— 0) is also a solution. 
Thus, having established the existence of one front (c_|_, 0+) with c+ > 0, we may take 
c_ = ^c_|_ and 4>~ = 4'^{—0- This situation, sometimes referred to as "gluing a 
front and back together" is typical for the kinds of scalar equations with comparison 
principles for which (HI) and (H2) have been established. 

We now mention some examples for which a comparison principle has recently been 
leveraged to obtain stability. Consider the following convolution model for phase tran- 
sitions 

U7i = ^ JkUn-k - Un + f{uk) 

with / bistable. Existence of traveling fronts was established in [5] under an ellipticity 
assumption on the convolution kernel J; asymptotic stability was established in |22j . 
Chen and collaborators studied 

Un = ^ an,kUn+k + /(Ufc) 
— n|</co 

in the case that the kernel an,k is periodic in n and elliptic, and the nonlinearity / is 
of bistable type [5]. In both of these examples, the authors do not verify (H2) directly. 
However, their results imply (H2). 

Another situation for which results exist is front-stacking. In any of the above 
examples, we can replace the bistable nonlinearity / with a tristable nonlinearity, e.g. 
g{u) = u{u + 1)(m — 1)(m — ai){u — 02) with —1 < ai < < 02 < 1. We can restrict 
attention to u e [—1,0] and apply the above results for the bistable case to establish 
the existence and stability of a monotone front (c_,(/)~) connecting —1 to 0. Similarly 
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we can restrict attention to m G [0,1] to obtain a second monotone front (c_|-,0+) 
connecting to 1. In this case Theorem 11.21 establishes the existence of a monotone 
solution connecting —1 to 1 with a long plateau at which grows longer over time. 

One situation to which our results do not apply is front stacking in conservation 
laws. This is because in conservation laws there is a line of equilibria at the constant 
solutions which generates an additional neutral eigenvalue, violating (H2). Stability for 
fronts in semi-discrete conservations laws was established in [1]. However, the presence 
of an additional neutral mode complicates the analysis both for the stability of a single 
wave and for the interaction; this lies beyond the scope of this paper. 

We now describe several models for which strong numerical and analytical evidence 
exists for (HI) and (H2). 

Vainchtien and Van Vleck [28] derived the model 

Vn = -^{zn + Zn~i - 2y„) - /(y„) 

Zn = -^{yn + Vn+l - 2z„) - /(z„) 

in the study of martensitic phase transitions. The state variables x and y denote even 
and odd lattice sites in a chain with both nearest- and next-nearest- neighbor coupling. 
In the case that f{u) is the McKean sawtooth caricature, (HI) was rigorously established 
in [28| . In the same paper, a combination of analysis and numerical experiments strongly 
suggest stability for the linear variational equation about the traveling wave, i.e. (H2). 
As a final example, we consider the discrete Fitzhugh-Nagumo equation 

{eiln = d{u„+i + Un-l - 2Un) + U„{Un - l)(li„ - fl) - U„ 
i'n = M« - bVn- 

A method for constructing pulses in the small e regime is described in [6]. In the same 
paper numerical experiments are reported which strongly suggest that these pulses are 
asymptotically stable. 

We remark, finally, that stability of pulses is generally more challenging than stability 
of monotone fronts in PDEs as well. This is because monotone tools such as the Krcin- 
Rutman theorem are not available to control the location of the discrete spectrum. 
Instead, the spectrum is usually controlled via Evans function methods e.g. [21]. The 
Evans function is built on top of exponential dichotomies for the spatial dynamical 
problem (e.g. ((B])) and requires finite-dimensional unstable manifolds. In the continuum 
case (O becomes an ODE and this is not a problem. Exponential dichotomies have been 
constructed for mixed type equations such as © (see e.g. [U]). However, the unstable 
manifolds are typically infinite dimensional. The appropriate generalization of Evans 
function techniques to equations like (O is an area of active research [3]. 

2 CMS-type decomposition to stable and center di- 
rections 

In the study of stability of traveling waves for PDEs it is standard to change coordinates 
to a moving frame in which the traveling wave becomes an equilibrium. Lattices do not 
admit such a moving frame. Nevertheless, traveling waves on lattices are shift-periodic. 
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that is (j){n — cTc) ~ (j){n — 1) when Tc = ^. The stabihty theory for travehng waves on 
lattices, developed in [2, is based on a Floquet theory for the time Tc = - map. A key 
step in the development of this Floquet theory is the construction of local coordinates 
which separate the neutral mode associated with translations of the traveling wave 
from the rest of the phase space. The purpose of this section is to develop a similar 
decomposition for the situation when two traveling waves are present. 
Let 



'-{■~c±t) e x^b 



and 



t e 



e CiR, GL{X^b)) with the following 



'0 - {P^it) 
Lemma 4.1 in [3] shows that there exist Z± 
properties, which hold for all 6* G M: 

• Z±{0) = Id 

. z±ie + i/c±)^sz±{9) 

. Z±(0)p±(O)=p±(0) 

Note that in [5] the authors work in spaces Z^, which have norms which are invariant 
under the shift S, and thus they can conclude (by the second property) that the operator 
norm of Z±{9) is bounded independent of 6. Our spaces Xb are not shift independent 
and thus the operator norm of Z±{6) may be large if 6 is large. See below. 

Now fix codimension one subspaces E'^ C Xz^b which do not contain p± and define 



X El 



by 



,y±)=p±(0±) + Z±(0±)y^ 



Proposition 4.2 in [9j ensures that {6^,y^) can used as local coordinates nearby V^'^, 
where the chart is given by 

Letting X±{t) = {9^ (t) , {t)) we now derive equations of motion for 6^ and 
y^. We carry the details out for the "— " component. Differentiating X± with respect 
to time and using ([3]) gives: 

z- {0- ) { [r (0) + q- {e- )y-] e- + y-}^F_{x-,x+) 

where the operator valued function is given by 

q-{e) := Z_{e)-^DZ_{6). 
Note that q^ {9 + —) = q {0) and thus the operator norm of q^ is bounded uniformly 



After multiplying both sides hy Z-{0 ) ^, apply the functional v & Xj^ , defined 
so as to annihilate Ej_ (and thus y~) and which maps p~{0) to one. This yields 

1 



= 9" 



,y 



We can solve for y 

r ^Y~{e-,y-,l 



:= , ^ , {Z^i9)-'F^{X-,X+)) . (7) 

1 + v [q [0 )y ) ^ 

) := z_{0-)-^F_{x-,x+)-[m + q-{o-)y-]e-{o-,y-,o+,y+) 



Id- 



p{O}+q{0-}y- 



l + iy-(q-(e-)y-) 



K-)J Z^{0-)-^F-{X-,X+). 

(8) 
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(10) 



Similarly, we choose E A'*^ which annihilates y+ and maps p^(0) to one to derive 
similar equations for and . 
Define 7±(i) ■.= e^{t)-t and 

r-{j-,y-,-,+ ,y+,t) := Q- {-f- + t,y- ,j+ + t,y+) - 1 

= il + ,.-{q{e-)y-))-' [u-Z^i9-)-'F^{X-,X+) 

- il + ,y-{q{9-)y-)),.-Z^i9-)-^F^ip-ie-),0)] 

(9) 

and similarly for r+. In the second line we have used the fact that 

1 = ly-ip-m = iy-{z(6-)-'r{e-)) = v-{z{o-)-^F^{p-{e-),o)). 

Therefore ^ becomes 

y- =Y-{e-,y-,e+,y+) y+ = ¥+{9- ,y- ,6+ ,y+) 
7" = r-(7",y",7+,?/+,i) 7+ = T+{j- ,y- ,'-f+ ,y+ ,t) 

Now we let 

Y_oij,y,t)^Y^it + j,y,0,0) 

and 

Y^ii-f-,y-,-f+,y+,t) = Y^{t + j-,y-,t + j+,y+)-Y^o{j-,y-,t) 
and similarly for Y+. Let A±{t) := DyY±a{t,0). Then (ITO| becomes 

y- = ^-(%- + {{A-{t + A-{t))y-} + {Y^oh-,y-,t) - DyY_oil~,0, t)y~} + ^-1(7^, y", 7+, y+, t) 

y+ = A+{t)y+ + {{A+{t + 7+) - A+{t))y+} + {Y+oh+,y+,t) - i?,y+o(7+, 0, t)y+} + Y+i{^- , y- , j+ , y+,t) 

7+ = r+(7-,?/-,7+,?/+,t) 

(11) 

This system is equivalent to ([31) in a neighborhood of Vq x Vq'" . 

3 Estimates for the right hand side 

In this section we prove a series of useful estimates for the right hand side for (|lip . The 
most important term is Z-{9~)~^[F{X~ , X'^) — F{X~,0)] which appears in both y_i 
and r_. 

Proposition 3.1. 

[{F{x-,x+) Fix-,0))] \u < ^^^^KlK^l (1 + I'y^il'^^^ + lly^ll^- + \\y^\\^-^\\y'\\^^) 

(12) 
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Proof. We compute 

Z-ie-)-^ [iF{X-,X+)-F{X-,0))] =Z-{0-)-^ [H- {G{X~ + X+) - GiX-) - G{X+))] . 

Note the following consequence of the mean value theorem. Recall that G{X) is of the 
form G{X)n = .g((7ViX)„, • • • (iVjX)„). Let x denote the vector {NiX+ , ■■■ , NjX+) 
and let y denote the vector {NiX" ^ • • • , NjX~). Use the mean value theorem to write 
g{x + y) — g{x) = Dg(x + ty)ydt and g(y) ~ Dg{ty)ydt so that, after using the 
fact that Dg is locally Lipschitz, we obtain 

\g{x + y) -g{x)-g{y)\ 

= 1 / {Dg(x + ty)-Dg{ty)}ydt\<C\x\\y\ 
Jo 

<GJ2\N^X~\\NkX+l (13) 

where the constant C may be chosen uniformly on bounded sets of x and y. 

Now let [tJ denote the greatest integer less than r. Notice that the second property 
of Z± implies that 

z±{e^) = z±{[e^c±\/c± + e^) = s^^^z^ie^) (i4) 

and 

where m± = [6*^0^] and 6 = c±9 — m± e [0, l/|c±|). Since 6 is restricted to lie in a 
compact set and 9 Z±{9) is continuous, it follows that there is a universal constant 
C such that the operator norm of Z±{9) and its inverse are bounded by G. 
This together with implies: 



\\Z^{9-)-' [{F{X-,X+)-F{X-,0))] \U. 



i,k 


\s- 


771 _ 


H- 


{\mx- 


-\\NuX+\]\U, 




i,k 


\\S' 




H' 


'{N^X- 


]\\S-"^-H-{N,X+}\\\^^ 




i,k 


\\s- 


-m 


H 


-{N^X- 


}\hA\S'"'-H-{NuX^}\l^ 




^^E 

i,k 


\\s- 


-m 


H 




'fi-M\xA\S''"'-H-{Nup+{9+)}\\l^ 


(15) 


+^E 

i,k 


\\s- 


-m 


H- 


-{N^p- 






+^E 

i,k 


\\S' 


-m 




'{N,Z- 


i(^-)v~}\hA\s^'''-H-{Nup+{e+)}\t^ 




+^E 

i,k 


\\S' 


-m 




-{N,Z^ 


i(^-)y'}\\\x,\p-'''-H'{NkZ+{e+)y+]\\\^^ 





In the second line of p5|) we have used the fact that H (xy) = {H x){H y). In the 
third line we have used the estimate ||C^^||a'6 < 11^^11 a'JI^IU"' • Each of the last four 
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terms corresponds to one of the four terms on the right hand side of the estimate in the 
proposition. 

We first estimate the contribution from the rightmost pulse, which is small be- 
cause of the cutoff function H^. 

|15-'"-i/-7VfcP+(0+)|U=o 
= sup(l — h{n — ct + m-))\Nk4''^ {n — c+9+ + m_)| 

nGZ 

<l|A^fc||£(;t_.)ll0+IU-. sup ((1 - h{n - ct + m_))(l + e-K«-c+('++™-))-iA 

<C(l + e-''(^*-'=+^+')-i 
<C(1 + e-''(-^^'^*-=+T+))-i 

pC+|7+ 

<c- 



1 e2('=+-'=-)* 

We now estimate the contribution from y"*". 

\\S-'''-H-NkZ+{e+)y+\\f^ 
= \\S-'^-H-NkS"'^Z+{e+)y+\U^ 

= sup(l - h{n -ct + ■m^))\[NkZ+{0+)y+]{n + m_ - m+)| 
nez 

<\\Z+ie+)\\ci;,_jN,\\c^;,_Jy+\\;,_, sup ((1 - hin - ct + m_))(l + 

<C||y+|U_,(l + e-''(^"*-"^))-i 
<C||y+IU_,(l + e-''(^"*-'=+^++^~+))-i 



e' 



1 + e2(=+^=-)* 



<c||y+IU_, 

We have used (|14p together with the fact that S commutes with N^- 

We now estimate the contribution from the leftmost pulse p" , which is bounded: 

\\S-'^-H-N,p-{9_)\\x, < !|7V,5-™-p-(0_)|U, < \m\cix,)U-i-c-l'+9-)\\x, < C(l+|7-|). 

Finally, the contribution from : 

\\S-^-H-{mZ-{e-)y-}\\x, < \\S-"'-N,S"'-Z^{0-)y-\\x, 

< \m\cixJ\Z-{e~)\\cixJy-\U 

< C\\y~\U. 

In the first line we have used the pointwise boimd |(i/^{X})„| < \Xn\ and in the second 
line we have used the fact that S commutes with Ni. Arranging the estimates for p^ 
and y^ completes the proof. □ 

Having estimated this crucial term we are now ready to bound the right hand sides 
of the evolution equations pT|) 
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Proposition 3.2. 

(16) 
(17) 

Proof. Wc first estimate 

|r_| <c{\\z^{e-)-'[F{x+,x-)-F{x+M\U + \\z^{0-)-'[f{x-,o) - ^^(p-,o)]||;,, + \\y-\u) 

Here we have used ([TT]) together with the fact that the terms iq^{e-)y-) ^^"^ 
\\q-iO^)\\c(^Xi-'i are bounded uniformly by a constant. We now estimate the term 

\\z^i0-)-^[F{x-,o)-F{p-,o)]\u = \\z^ie-)s-"^-[F{x-,o)-Fip~,o)]\\ 

= II Z_ {§- ) -1 [F(5-™- (p- + 5'"- Z_ (9- )y- , 0) - F{S-^-p- , 0)] |U, 

< cily-iu,. 



Combining this with Proposition 13.11 yields (|T6|) 
We compute 



l + v-{q-{e-)y-) 



Thus 



l>^-ilU. < ||i-f+P^^(^z.-(.)||,^||z_(^?-)-MW+,x-)-F(xH 

C(l + l7-|)el-+-'+l) ^1 , |L,-|| | 



Here we have used the fact that the operator norm of q- is bounded uniformly in 6 and 
that ||y~|| can be made small to bound the first term and Proposition 13 . 1 1 to bound the 
second term. 
Since 

||(^_(t + 7-) - A-{t))y~\U < C|7-|lb-|U. 
and F_o(7^, 0, t) = 0, we also have the estimate 

\\Y^oil^,y^,t) - DyYM^^,0,t)y-\\x, < C\\y-\\%^. 

In light of (fTT|) . this yields p7|) and hence completes the proof. □ 
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4 Proof of Theorem 11.2 



Proof. Let B±{t, to) denote the evolution operator associated to y — A±{t)y. It follows 
from (H3) and statement 1 in Theorem 5.3 in '^W that 

for some C > 0, A > 0, independent of b. 

After applying the Duhamel formula to the equations for y and using Proposition 
13.21 the equations for y^ in pT|) give: 



lly-WL^ < ||B-(t,to)||£(;t,)||y"(^o)|U^+ / \\B^it,s)\\^^^^^\\y-i.s)-A{s)y~{s)\Uds 

J to 

Jta 

]^ -|- g 2 ^ + <-- 

(18) 

Similarly, for 7" we have: 

t 



b^Wl <C / (||y-(s)||^^ + |7-(s)|||y-(s)|U, 

+ ^ \ , , ^ (1 + \\V-{^)U. + + ll2/+(^)lk-J|jy'(^)lk.) rf^- 

(19) 

There are similar estimate for ?/+ and 7+. 

Now let 5 := ||y-(to)IU, + b+(io)IU,' let := |(c+ - c_), let a := min{A/4,6*} 
and define 

Kt:^ sup f|7"(i)| + |7+(i)|+e°(*"*"H\^+e-"")-'(||2/-(i)L •) 

to<t<T L' V" / 



to<t<T 

Thus 



whenever < t <T . Note that = and Kx is increasing with T. (This increase is 
continuous since our LDE is locally well posed.) Our theorem is proven if we can show 
that Kt is bounded uniformly for all T > Iq. We choose T so that Kt < 1- 
Then iW) gives for to < ^ < 7": 



y~{t)\\^^<CSe-^'^'~''>^ +CK^ f e-^^'-'^^e-^'ds + C I e'^^*"'^) ^ 

^ ./ J to 1 + e 



-ds. 



for some constant C which is independent of T, 5, and to- Here we have used the 
fact that one dominates ||2/^(t)|| a"^;, and that + *" dominates (v^ + e~'' 
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Integrating the exponentials and using the fact that for s large enought, j^_^_^2b* s is well 
approximated by e"^** * we obtain 

ga(t-to)(^+gb-to)-l||y-(i)||^^ < C e-^^-'^)^*-*") + A'2 + (g-(2b*-a)(t-to) + g- ( A-a) (t-to) ) 

We can control the right hand side of in much the same fashion, though we 
omit the details. All together we can show there exists C* > (independent of T, S, 
and to) so that: 

Kt <C* {y~5 + K^ + e"''**") . (20) 

There exists positive constants tg, and < K* < 1 so that if < ii' < K* , Q < 5 < 5q 
and t > to, then 

C* ((5o + K^ + e^****") < K/2. (21) 

Let r* be the smallest time greater than to for which Kt = K*, if such a T exists. 
Otherwise set T* = +oo. Notice that if T* ~ +oo that we are done with the proof 
of Theorem 11.21 {to in this formulation corresponds to ^^'^^ in the statement of the 
theorem) Suppose that T* < +oo. If so, then (^0)1 and (PT|) imply that 

Kt < Kt/2 

which is a contradiction. We are done. □ 
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